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Abstract 

In this article we focus on estimating the quadratic covariation of continuous semimartingales from discrete 
observations that take place at asynchronous observation times. The Hayashi-Yoshida estimator serves as 
synchronized realized covolatility for that we give our own distinct illustration based on an iterative syn- 
chronization algorithm. We consider high-frequency asymptotics and prove a feasible stable central limit 
theorem. The characteristics of non-synchronous observation schemes affecting the asymptotic variance 
are captured by a notion of asymptotic covariations of times. These are precisely illuminated and explicitly 
deduced for the important case of independent time-homogeneous Poisson sampling. 
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1. Introduction 

Nonparametric estimation methods for the quadratic variation of semimartingales have become an issue 
of great interest in recent years. One reason is the interpretation of the quadratic variation of the continuous 
part as integrated volatility in financial modeling. 

If a semimartingale is observed discretely at times t^, < « < n on a finite time horizon [0, T], the sum 
of squared returns (increments of the semimartingale), called realized volatility, converges to the quadratic 
variation as sup {ti — ti-i) — as n — oo. The same fact pertains to the multi-dimensional case where 
the realized covolatilities of two processes converge to the quadratic covariations. More usually multivari- 
ate data, in particular financial time series, are recorded at times following non-synchronous observation 
schemes. Therefore, realized covolatility estimates most commonly incorporate a previous-tick interpo- 
lation approach. Though, this machinery leads to the so-called Epps effect |8| that realized covolatilities 
tend to zero as the sampling frequency increases. Especially for the more and more available ultra high- 
frequency financial tick-data this issue poses problems. 

A solution for the asynchronous estimation problem has been proposed in ifTSl . We call this estimator 
which arises as realized covolatility from all products of returns with overlapping observation time instants 
Hayashi-Yoshida estimator. Our investigation of that estimation approach leads to several useful rewrit- 
ings and interpretations. The final representation is based on an iterative synchronization procedure which 
has been used first in |19|. This synchronized realized covolatility and the data aggregation technique for 
synchronization can serve as a basis for combined approaches in various generalizations of the underlying 
statistical model. A very important enhancement of the model in that we take market microstructure noise 
into account is covered in |01 and |[5| by extending the synchronized realized covolatility to a generalized 
multiscale estimator. 

The asymptotic theory developed in the article on hand is grounded on stable limit theorems for semi- 
martingales from lfT6l . We obtain a stable limit theorem for the process associated with the estimation 



'Financial support from the Deutsclie Forscliungsgemeinschaft via SFB 649 'Okonomisclies Risiko'. Humboldt-Universitat zu 
Berlin, is gratefully acknowledged. 



1 



error of a Hayashi-Yoshida estimator for the quadratic covariation at time t E [0, T] and for the overall 
estimator by the marginal distribution at t = T the stable weak convergence to a centred mixed normal 
limiting distribution. The random asymptotic variance splits up in two terms induced by an idealized syn- 
chronous approximation and an additional error due to the lack of synchronicity. 

The article is arranged in six sections. In the following Section[2]we give insight into the concept of stable 
weak convergence and a short review on the essential theory from the literature. In SectionjSlthe Hayashi- 
Yoshida estimator and our related synchronization algorithm that we have first presented inl4l is revisited 
and the asymptotic theory including the key result is provided in Section]?] The detailed proof of the central 
Theorem ]2] is postponed to the Appendix A Following some simple illustrative and motivating examples 
before, Section ]5] comes up with the analysis for the important time-homogeneous independent Poisson 
sampling case for which we evaluate all ingredients of the asymptotic variance explicitly. To benefit from 
the stable central limit theorem and provide a basis for statistical inference, we give a consistent estimator 
for the asymptotic variance in Section ]6] 



2. Stable convergence and Jacod's stable limit theorem revisited 

This section is devoted to the notion of stable weak convergence which will be an essential concept 
for the development of our limit theory throughout this article. The concept of stable convergence goes 
back to ||2TI and results about stable limit theorems were extended in HI and |l9|. The reason what makes 
stable weak convergence a key element of our asymptotic considerations, is that it allows to conclude joint 
weak convergence when we derive results about asymptotic mixed normality. In the case that a sequence 
of random variables weakly converges to a mixed Gaussian limiting random variable VZ, with Z 
being standard normally distributed, Z ^ N(0, 1), and a strictly positive random variable V, independent 
of Z, we cannot derive confidence intervals if the distribution of V is unknown. However, if a consistent 
estimator for the asymptotic variance is available (in the sense that — ^ V^), the stable weak 
convergence will assure that (X„, V^) {VZ, V"^) jointly and also that XnjVn ^ Z. The last implica- 

si 

tion also holds in a stable version: XnjVu Z. For that reason, if we are in situations as described above 
we gain from proving stable weak convergence which paves the way towards statistical inference. 
Next, we present the formal definition and the main properties of stable weak convergence of sequences of 
random variables. 



Let (Xn) be a sequence of random variables defined on some probability space (fi, T , P) and taking values 
in a Polish space (E, £). We say that the sequence converges weakly in to X if for any bounded 
random variable Z 

lim E [ZX^] = E [ZX] 

n— ^oo 

holds. 



Definition 1. For a sub-cr-field Q C J^the sequence of random variables is said to converge ^-stably, 
if there is a random probability measure ^ovi{Q, y. E,Q ® £) such that 



lim E[Z/(X„)] = / fi{du,dx)Z{uj)f{x) 
JnxE 

for all / G Cb{E) (continuous and bounded) and tj-measurable bounded random variables Z. 

st 

If G ^ J^, we say (Xn) converges stably in law to X (X„ X). 

Remark 1. t^-stable convergence is the weak convergence in of E [/(X„)|C7] for all / £ C{E) to /io /. 
This implies convergence in distribution to the probability measure i' defined by 



v{B) = J ^i{dio,B)l 

{X (uijeBj^idi^) 

If (Xn) converges stably, the limiting law is fj.{ft, ■ ). 



2 



The following proposition states some useful equivalent characterizations of (t/ — )stable convergence. 
Proposition 2.1. (X„) converges Q-stably is equivalent to: 

(i) For every Q-measurable random variable Z on Vt, (Z, X„) converges in law. 

(ii) For every Q-measurable random variable Z on f2, {Z, X^) converges Q-stably. 

(iii) The sequence is tight, and for allG (z Q and f G C{E), the sequence E [1(5/(X„)] converges. 

This proposition is proved in ifTTl as part of Proposition IX. 1.4. Stable (weak) convergence is a stronger 
mode of ordinary convergence in distribution. It is weaker than convergence in probability, but we empha- 
size that the limit depends on the limiting random variable X itself and not only on the distribution of X. 
If [Xn) converges stably to X, X is defined on an extension (51', F\V = fj.) of the original probability 
space, so that 

V/ e CiE) : lim E [/(X„)Z] = E' [f{X)Z] . 

In the situation that we face in this article, a Gaussian random variable which is independent of F will 
appear as limiting law. In this case we call the extension of the original probability space orthogonal. In 
particular, if Xn ^ X holds with an J^-measurable random variable X (L^— convergence), the foregoing 
proposition yields that (X„, X) {X, X) and, hence (X„ — X) holds, which implies convergence in 
probability. Therefore, in all cases where stable weak convergence is a suitable adequate concept, limiting 
laws are defined on a genuine extension of the original probability space. 

The following proposition gives the result that stable convergence is the suitable concept to derive feasible 
central limit theorems and confidence intervals if the asymptotic variances in limit theorems are unknown 
random, but can be estimated consistently. 

st 

Proposition 2.2. Let {Xn, Vn) be real-valued random variables defined on (fl, P). If Xn X with a 
mixed normal limiting random variable X ^ N(0, V'^) and Vn — ^ V with V being T-measurable. Then 

Xn/Vn ^ N(0, 1) 

holds true. 

Note that we use the same denotation expression for mixed normal laws and common normal laws 
and the difference becomes clear out of the context and by the specific variances. On the assumptions of 

st 

the proposition (X„, Vn) {X, V) is implied and the convergence of Xn/Vn follows by the continuous 
mapping theorem. This proposition is part of Proposition 2.5 in ll20l . We restricted ourselves to real-valued 
random variables in the last proposition. A more general version can be found in 1 17 |. 
The concept of stable convergence also carries over to stochastic processes. For this extension of stable 
convergence to stochastic processes, or more precisely to semimartingales, the Polish space E in Definition 
[T|is chosen to be the Skorohod space. The following limit theorem for stable convergence of continuous 
local martingales will be the foundation for our later deduced limit theorem in this article: 

Theorem 1 (Jacod's theorem: A martingale version). If{Mt,J-t) with < t < oo is a continuous local 
martingale defined on the probability space {ft, T , P), we denote by M.^ the set of bounded {Tt)-adapted 
martingales orthogonal to M = {Mt,J^t) what means that [M, M^] = 0. If {X"') is a sequence of 
continuous {Tt)-odapted local martingales for which 

[x\~^Vt yt e [o,oo) (1) 

with a continuous process V holds, the following two conditions 

[X",M]j -A Vi e [0,c3o) (2a) 

[X",iV]j ^0 yt e [0,oo) andyX e (2b) 

are sufficient that (X") converges {T)- stably in law to Wvt, where W is a standard Brownian motion 
independent of T. 
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This theorem is a simplified martingale version of the more general theorem 2-1 in lfT6l . A similar 
special version of the theorem is also used in fTOl. A comprehensive illustrative overview on Jacod's stable 
limit theory and further motivation and applications of this result can be found in li20J . A discrete-time 
version of that theorem (cf. 3-1 in fl6l ) is the following: 

Corollary 2.3. Assw7ie that Z"^ = .<(^n,i is the endpoint of a discrete martingale and the Xn^i 

are Tt„ ^-measurable square integrable random variables and {Wt^ Tt) a Brownian motion and ATn^i — 
Tn,i+i — Tn.i — >■ fli n — >■ oo. If there exists a predictable process {vs)s>o such that 



E 

T„ i<t 



EK.|^t„,._J ^ / vlds, (3a) 
J Jo 

Ve > : ^ E [Xl^ l{x„,.>e} I-^t„_ J ^ , (3b) 



.,i<t 



E [X„,,(Wt„,. - VFt^^.^IJ-t^^.J ^ , (3c) 



E 



E [X„.(Mt„„ - Mt„_J|J-t„„_J ^ , (3d) 



for all bounded Tt-martingales with Mq — and [W, M] = 0. Then the following stable convergence of 
the process holds true: 



Z'^^Zt= f VsdW^ (4) 





where is a Brownian motion defined on an orthogonal extension of the original probability space. 

The limiting process in the foregoing Theorem [T] is a time-changed Brownian motion. The Brownian 
motion is of central importance in the theory of continuous local martingales, since every continuous local 
martingale Mt corresponds to a Dambis, Dubins-Schwarz time-changed Brownian motion ■ For each 
(X") we have a representation as Dambis, Dubins-Schwarz Brownian motion W^„] and the sequence 
converges weakly to a limiting Brownian motion Wy by the asymptotic Knight-theorem. We refer to 
Theorem 7.7 in [6| for a proof The conditions ( |2a| i and ( |2b] i about the quadratic covariations converging 
to zero in probability ensure that the weak convergence to Wy is stable. 

For one fixed < T < oo we have the result that XJ^ converges stably in law to a centred mixed normal 
distribution: 

X^^N (0, Vt) ■ (5) 

The independence of the limiting Brownian motion W and (V, Y) for any J^-measurable random variable 
Y assures that {Wvt,Y) has the same law as (VtZ, Y) with Z N(0, 1) and independent of {Vt, Y). 
Note, that in the original theorem 2-1 in lfT6l for semimartingales the same conditions as in our Theo- 
rem [T] are imposed for the predictable quadratic (co-)variation processes that coincide with the quadratic 
(co-)variations for continuous semimartingales. Additionally, a condition that the drift can be neglected 
asymptotically is imposed. Compared to Theorem 3-1 in lfT6l . we allow for non-equidistant discrete par- 
titions which does not harm the deduction of Theorem 3-1 from Theorem 2-1 in [16|. A conditional 
Lindeberg-condition ( (3b] i and a convergence condition on the conditional variances (|3a]) are analogous as 
in central limit theorems for triangular martingale arrays. The main difference to the stable limit theorem 
Corollary 3. 1 in fT2l (page 58 ff. ) is that a certain nesting condition on the filtrations is replaced by con- 
ditions ( (3c] i and (3d|. Usually the reference Brownian motion W is given and "fully generates" the X„ jS 
in the sense that (3d]) holds. 
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The theorem also extends to a muhi-dimensional setting which is formulated separately in the next corol- 
lary. For this purpose let M* denote the transpose of a vector M and the {d x r)-dimensional quadratic 
covariation [AI,N*]^ ■ with 1 < i < o? and 1 < j < r for a d-dimensional AI and 

r-dimensional N. Recall that convergence in probability of a vector is equivalent to convergence in proba- 
bility for every component. 

Corollary 2.4. Let {Mt,J-t) be a d-dimensional continuous local martingale and M.'^ again the set of 
{J-t)-adapted bounded martingales orthogonal to M (to all components). A sequence of r-dimensional 
continuous {J^t)-adapted local martingales (^") with 



where is a predictable JV (g) R'' process, and 



Vt 



WgW* ds 



(6) 



(7a) 



Vi G [0,oo) and^N e 



(7b) 



converges stably in law to the process WsdWs, where W is a r-dimensional standard Brownian motion 
independent of J-. 

Jacod's theorem provides a convenient stable central limit theorem for our purpose. Nesting conditions 
on the sequence of filtrations that are required for other stable limit theorems as in 1 12 | and |22| are not 
satisfied here. 

Furthermore, the concept of stable convergence enables us to prove the stable weak convergence to mixed 
Gaussian limiting random variables under an equivalent martingale measure P after a Girsanov transfor- 
mation, where the drift processes are zero. Stable convergence guarantees that the asymptotic law carries 
over to the case with drift under the original measure P. It is in this sense commutative with measure 
change (cf. fTSl). If we have the result that Z„ ^ to + AVAR • N(0, 1) under P with a standard Gaussian 
distribution independent of T, defined on an orthogonal extension of the original probability space and 
J^-measurable bounded random variables m and AVAR, the same convergence holds true under P. Since 

st 

stable convergence Z„ Z implies for all / e C(£) and J^-measurable bounded random variables X 

E [X/(Z„)] - E [(dP/dP)X/(Z„)] ^ E' [(dP/dP)X/(Z)] = E' {Xf{Z)\ , 
by uniform integrability of X/(Z„)(dP/(ilP) with 



dP/dP = exp 



^ -isdB, + -J^ ds ) , 



where ctsTs + Ms = 0. 



3. A synchronized realized covolatility estimator 



Assumption 1. On a filtered probability space {fl,J-, [Ft) ^V), X — (Xt)teK,+ cind Y 
continuous semimartingales defined by the following stochastic differential equations: 



{Yt)t(iK+ are 



dXt 
dYt 



pf dt + dB^ 
pj dt + o-f dBY , 



with two (Tt)— adapted standard Brownian motions B^ and B^ and ptdt = d \B-^ ^ ^'^^ drift 

processes p^ and p( are {J-t)-adapted locally bounded stochastic processes and the spot volatilities 
ij^ and aj and pt are assumed to be {J-t)-adapted with continuous paths. We assume strictly positive 



volatilities and the Novikov condition E 



< oo for X and Y. 
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We consider the estimation of the quadratic covariation [X, F]y of two continuous semimartingales, 
also called Ito processes, X and Y as defined in As sumption [T] from discrete observations following non- 
synchronous sampling schemes. 

We impose the following regularity assumptions on the underlying asynchronous sampling schemes: 

Assumption 2. The deterministic observation times 7"^'" = {0 < tg"^ < t^"^ < . . . < ti"^ < T} of X 
and 7"^'™ = {0 < Tq™-* < t[^^^ < . . . < r/,™'' < T} ofY are assumed to be regular in the following 
sense: There exists a constant < a < 1/3 such that 

6^^ sup ((4-^-t^Xt^-\T-t^-A =0(n-V3-a)^ (8^) 
ie{i,...,7i} vv / / \ / 



(8b) 



We consider asymptotics where the number of observations of X and Y are assumed to be of the same 
asymptotic order n = 0{m) and m = 0{n) and express that shortly by n ^ m. 

For synchronous data n — m and — Tj*^"^ for all i e {0, . . . , n} holds. In the non- synchronous case 
the number of observations (n + 1) of X and {rn + 1) of F may differ and the sets of observation times 
^x,n ^Yso contain times ^ ^Y,m ^j^^ ^j") ^ q'^^^^. We work within the general model where also 
synchronous observation times can take place and hence 7~^>™ and 7"^ " are not assumed to be disjoint. 
In the following, we omit the superscripts (n) and (m) for observation times to increase the readability. 
Although the sequences of observation times are modeled deterministically, we remark that the case of 
random sampling times that are independent of the observed processes is included in that analysis regarding 
the conditional law given the observation times. 

We use the short notation AXt . ,i = 1, . . . ,n from now on for increments Xt^ — Xti_^ and analogously 
for Y . In ifTsl the consistency of the estimator 

(HY) 

[X,Y]rj. = 2^2^ AXt^AyT-^l[niin(t.^.rj)>max(t,_i,T3_i)] , 

i=l j = l 

is proved, where the product terms include all increments of the processes with overlapping observation 
time intervals, for a similar model of discretely observed Ito diffusions with deterministic correlation, drift 
and volatility functions. Consistency directly carries over to our setting including random correlation, drift 
and volatility processes. The estimator is also in our setting, furthermore, unbiased if drift terms are zero 
and else asymptotically unbiased. In [14J it has further been shown that on stronger regularity assumptions 
on the observation schemes this Hayashi-Yoshida estimator is asymptotically distributed according to a 
Gaussian law. 

For a general strategy leading to a synchronization mechanism that keeps to the Hayashi-Yoshida approach 
and its valuable properties, we focus on an alternative useful method to handle the asynchronicity of the 
data. It has been introduced in llT9l . where it was called pseudo-aggregation. The method translates the 
Hayashi-Yoshida estimator into an iterative algorithm that allows to rewrite the estimator without indicator 
functions. This can be done by aggregation of addends for which partial sums are telescoping. A first 
simple rewriting of the Hayashi-Yoshida estimator is obtained by taking the sum of the products of aU 
increments of X with the telescoping sums of aggregated observed increments of Y for that observation 
time instants overlap with the according observation time instant of X (or in the symmetric way): 

(HY) " ^ / ^ 

[X,Y]rj. =^AXfJ ^ AF^^l[niin(t.^^^.)>max(ti_i,T,_i)] 
1=1 \j6{l,...,Tn} 

m / 

= ^Ar^^ E ^'^t.l[nnn(t,,r,)>max(t._i,r,_i)] 

j=l \ie{l,...,n} 
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first step: 

• for to < To and /io = min (w e {1, . . . , n}\TQ < t^): 

H° = {to,...,VJ and g" = {To} 



if To < 

• for to = To: 

= {to} and - {to} 
gi = 1 and ri — 1 

• for to > To and wq = min (/ e {1, . . . , m}|to < t;): 

H° = {to} and g° = {To,...,T^J 



fwo + 1 if to = T, 

gi = 1 and ri = < -e 
^ \wq if to < T, 



Wq 
Wq 



ith step (given and C?* ) 



• for tg; < Tr; and fj,i = min (w e {gi + 1, . . . , n}|Tr; < tm): 

H*-{t,.,...,t^J and e' = {T,J 

'^^ ' = if T,. < t^, ^'+1 - + 1 

• for tq. = Tr- : 

W = {t,J and = {t,J 
qi+i = Qi + l and r^+i = + 1 

• for tg. > Tr- and Wi = min (/ e {r^ + 1, . . . , m}|tg. < t;): 

n'^{tg,} and a' = {t,.,...,t„J 

, 1 J I = + 1 if to = Tw 

= + ^ ^ ^: if t: < 

Algorithm 1: Iterative algorithm for construction of the joint grid from asynchronous data. 



Defining the next-tick interpolation t, + := mino<j<„i (tj|tj > tj) and the previous-tick interpolation 
ti._ := maxo<j<m (tj |tj < ti), the last expression can be illustrated 

2=1 

The algorithm which we will use is a more enhanced method to aggregate the data in an adequate way. 
For this purpose {N + 1) sets W and CJ' are constructed, where TV < min (n, to), each set including one 
or more than one observation time of X and Y, respectively. This method to construct a joint grid for the 
observations of the two processes is described by Algorithm[T] 

The Algorithm[T|that we have first presented in 14J stops after (N + 1) steps when the last observation time 
is reached. We pass over from the original observations to the sums of observed increments over sets 
W and over sets C/% respectively. The observations are grouped together so that the resulting realized 
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covolatility estimator 



A'' n ni 

i=0 1=1 j = l 

calculated from the 'synchronized' observations 

for the integrated covolatility wiU coincide with the one by ITSl stated above. We use a different illustration 
of this estimator compared to |fT9l making use of telescoping sums. 
With the denotation expressions from Algorithm[T] 

lii — max (fc|t/c e H*), Wi — max (fc|Tfe € Q^) and 

q, = min {k\tk eW), = min (fc|Tfe G g') ,ie{0,...,N} 

and for the purpose of a simpler notation 

Xg, = ^t,, , = ^'r^, , i € {0, . . . , iV} and 

with ?ci := Ao := tq, and F^' can be written as telescoping sums X^ = {Xg. — Xi.), Y^'^ = 
(Y^^ -Yx^) .This leads to 

[x,Y]^ = ^ (X,, - X J (y^, - n. ) , (9) 

4 = 1 

where summation starts with i = Ooii = 1 since the addend for i = is always zero. Although we use this 
specific new illustration throughout this article, we will call this realized covolatility of our synchronized 
observations also Hayashi-Yoshida estimator in the following. In this notation gi denotes the greatest and 
li the last observation time before the least element of the set W and analogously 7^ and A, of Q^. 



Example 

An illustration of the application of Algorithm [T] to observations is given in Figure [T] In this example, 
wehaveT^o = {to},g° = {to},H^ = {^i, ^2, ts}, = {ti},H' = {^3},^' = {t2,T3},H^ = 

{^4,^5,^6},^' = {Ti},n^ = {t6,tj},g^ = {rs},^^ = {t7,ts},g' = {Te},n^ = {is}, a' = 

{T,,Ts},n' = {tg},^^ = {t8,T9},H« = {tio},^' = {T9,rio} . 

The example highlights the important features of the synchronization procedure. The sets W and G'' are 
in general not disjoint and the maxima of consecutive sets can be the same time points. The minimum of 
a successive set can as well equal the maximum of the prevenient. Contrarily, consecutive minima are not 
equal. For further examples we refer to |,19J . Of course the example is just for illustration and the number 
of observations is much smaller than in practice. The synchronization of n + 1 = 11 and m + 1 = 11 
observations leads to iV + 1 = 9 synchronized observations in this example. 

The fact that we obtain {N + 1) < min (n, m) + 1 synchronized observations indicates heuristically that 
the efficiency of such techniques of covariance estimation mainly depends on the number of observations 
available for the less liquid process which is observed at a lower frequency. By Assumption |2] we restrict 
us to the case that n and m are of the same order. Thus for the suprema of times between two observations 

(iV-'/^'-") and 6^ =0 (iV-'/^"") 

holds with a constant < a < 1/3. 

In the next section, we show that on Assumption[l]and|2]the estimator (|9]l is V A^-consistent and, on further 
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Figure 1 : Example for synchi'onization using Algorithm[T] 

assumptions on the asymptotic behavior of the asynchronous sampling schemes, asymptotically normally 
distributed. Using standard interpolation methods such an estimator cannot be obtained. 
Another recent approach to deal with non-synchronous discrete observations in a general setting including 
market microstructure noise has been proposed by f2\. This method is also related to our approach. The 
so-called refresh times are the cumulative sums of waiting times until both processes are observed. As- 
sume that in the ith step of Algorithm[T]tg. < r^. holds. Then the next observation times of X are grouped 
together ending with the first observation time t^^^i < r,.. < t^. greater or equal than t^^. Then we start 
the next comparison step and compare this last observation time grouped to the set W to r^.+i, except 
for the case where two synchronous observations appeared, where we compare the two following times. 
Since in the completely asynchronous case at the refresh times only one of the two processes is observed, 
the refresh time method used in |2| includes a previous-tick interpolation for the unobserved process at the 
refresh times. Refresh times provide the 'closest synchronous approximation' to the asynchronous sam- 
pling schemes that we define in Proposition [TT|below. The number of refresh times which are denoted in 
this work hy Ti,i — 0, . . . , TV, equals the number of sets constructed by pseudo-aggregation. In a setting 
that also takes microstructure noise into account, a consistent estimator requires smoothing techniques to 
reduce the noise perturbation and the optimal convergence rate is slower (cf. |4|). The previous-tick in- 
terpolation, however, causes a negative bias due to asynchronicity when calculating the simple realized 
covolatility estimator based on the refresh time and previous-tick approach and it does not equal the esti- 
mator of Hayashi-Yoshida. The reason for this bias is that, due to the previous-tick interpolation, products 
of increments with overlapping observation time instants fall out of the realized covolatility. The pseudo- 
aggregation Algorithm [T] used in this work corresponds to the refresh time method when replacing the 
previous-tick interpolation by a next-tick interpolation for the right end points of refresh time instants. 
Then, the resulting realized covolatility of 'synchronized observations' 

(HY) 

i=l 
N 

= E (^T^ - ) {Yt-^ - (10) 

1=1 

coincides with the Hayashi-Yoshida estimator and has no bias due to asynchronicity. As figured out in 
the simulation study of |4| the asymptotically vanishing influence of the bias due to pure previous-tick 
interpolation also shows up in the setting with noise for finite sample sizes and mild noise variances for 
that combined estimators are constructed in ^ and f4], among others. 

Figure |2] visualizes refresh times T^, i = 0, . . . , 8 for our above given example. For this example the 
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Figure 2: Example for synchronization using Algorithm[T]including refresh times. 

realized covolatility calculated with refresh time previous-tick interpolated values equals 

(X,^ - X,J{Yr, ~Y,J + {X,, - X,,){Y,, - + {Xt, - X,,){Y,, - nj + 

{Xtg - XtJ(lVg - J + {Xt^^ - Xtg){YT^g - y^g) 

and is biased downwards due to non-synchronicity, whereas (|9]) yields 

(Xt3 - Xt,){Yr, - r.j + (Xt3 - Xt.,)(Yr., -Yr,) + {x^, - Xt,){Yr, - 
(X,, - x,j(y,, - i;j + - - + (x,, - x,,){y,, - 

which is an unbiased estimator for observations of processes according to Assumption[r| when drift terms 
are assumed to be zero. 

4. Asymptotic distribution theory 

In this section the elements for an analysis of the asymptotic properties of the estimator (|9]) are devel- 
oped where the emphasis is on the asymptotic distribution of the estimator 

The following technical Proposition constitutes the theoretical justification that the refresh times r> , 1 < 
i < N introduced in the foregoing section can serve as a convenient basis to decompose the overall esti- 
mation error of the synchronized realized covolatility (|9|. For every N these times induce a partition of the 
time horizon [0, T] that we call the closest synchronous approximation. 

Proposition 4.1. If we define T;^^' := min (5^, 7,), « = 0, . . . , A^, the set r^^"'^ {T^^^ T^^^^ 

induces a partition of the time span [0, T] in the sense that ij^[T^^\ T^'^^) = [Tq^\ T — T^'^). 
The following equality holds true: 

rr - min (,f \ 7f ') = ma. {l^^lAT^) ,^-l,.■■,N-l (11) 

and on Assumption^^dN :— sup,jg|]^ (^^^i^'' ~ ^i^i ) — ^ (A^^^/^^") holds. 

In the following we frequently leave out superscripts indicating dependence on A^ to guarantee clarity 
and increase the readability. 

Proof. Assume without loss of generality gi < 7^ for an arbitrarily fixed i e {1, . . . , A^ — 1}. Taking 
Algorithm [T] into account, we proof that ( fTTj l holds true. 

If gi < 7i, then the observation times 7,; and + :— min (ij. G T^\tk > gi) are compared in the (i + 
l)th step of the synchronization Algorithm [ll and gi^ = min (t^ G T'^ltk G 'H*^^) holds true. Thus, 
gi = li^i and ( [TT| i holds true. We remark that in this case 7^ e Q^^^ and thus 7^ > A^+i = 7i._ := 
max (rfc e T^\Tk < ji) > 7^-1. 

If gi = 7i, then the observation times 5^ and 7; + are compared in the {i + l)th step of Algorithm [T| and 
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= Ai+i = gi = 7i what implies ( fTT) . 
Equation ( [TT| does not hold true for i — 0,N and To = to A tq because we have set Iq = and Ao = To. 
Although consecutive maxima gi of the sets T-C and 7^ of the sets respectively, can be equal, Ti > T^-i 
holds for all i e {!,•■• ,-/V} because g^+i = (/i implies that 7^+1 > 7^ and 7^+1 = 7^ implies that 

> gi. Hence, the set 'T^v'^ induces a partition of the time span [0, T]. □ 

The times Ti,i = 0, . . . , N defined through ( [TT] i equal the refresh times from (2] as has been mentioned 
in the last section. We use Proposition 4. 1 to split the error of the estimator (|9]l for the integrated covolatility 



[X, Y]rp in two asymptotically uncorrelated parts. The error of the estimator (|9]l can be written 

N T 



where 



- J2 [i^T, - Xt,_,) {Yt, - Yt,_,) - ' PtafaY d^j (12) 
- / Ptcr^f^Ydt- / ptafaj dt 

Jo Jtr^AT,„ 

is a synchronous-type discretization error of the realized covolatility estimator evaluated with synchronous 
observations at the times Ti,i — 0, . . . ,N, which is the closest synchronous approximation to the asyn- 
chronous sampling scheme, and 



N 



1=1 

N 

+ J2 (Xt, - Xu) (K,. - ItJ 1{t,=<,j + {Xt^ - Xt^_,) {Yt^_, - l{T._i=i,} (13) 

i=l 

is the remaining additional error due to the lack of synchronicity. When we write the increments involved 
in the estimator (|9]) in the way 



{X,, - Xi^ ) = x+ + X/ + X- , (y,^. - n^. )^y+ + y;^ + y^ 



where X^ = Xg. — X^^ denotes the next-tick interpolation error at right-end points, X'^ — Xtj _i ^ Xi. 
the previous-tick interpolation error at left-end points, X^ = Xt^ — XTj_i , j — 1, . . . , N the increment 
over the time instant of the closest synchronous approximation and analogously for Y, and can be 
expressed: 



N 



1=1 

N 

= E i^^'^^' + + y^^ix^ + X-) + x-y/ + Y-xf) 



is an usual synchronous -type realized covolatility but incorporates an idealized sampling design at the 
times of the closest synchronous approximation for which we do not have observations in an asynchronous 
setting. Nevertheless, this idealized approximation turns out to be helpful for our further analysis. The 
error due to non-synchronicity A^ hinges on the interpolations that have to be carried out since we do not 
observe X and Y at the times Tf^ , 1 < i < N. The term is asymptotically centred since only products 
of increments over disjoint time instants remain whereas is an unbiased estimator for [X, Y]j,. Since 
either X or y is observed at a certain Ti,l < i < N, one of each interpolation errors in the illustration 
above equals zero. 
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to t, k h t, k k h k k 



Figure 3: Illustration of the synchi'onous approximation for our example. 

Proposition 4.2. The Brownian parts of and are uncorrelated. This means, that if we assume 
the drift terms to be identically zero in Assumption^^ and are uncorrelated. If the drift terms are 
non-zero, A^ and are asymptotically uncorrelated. 

Proof. A^ and are both centred. If Assumption [T] holds with /i^ = /i^ = 0, the expectation of the 
product of A^ and is zero, since the previous- and next-tick interpolated increments in ^T3[ are centred 
and uncorrelated to the other three factors in each addend of the inner sums. 

If we allow for non-zero drift terms, Assumption[T]and Assumption[2]ensure that the increments over time 
intervals due to the drift induce terms at most of order in probability by products of drift terms and at 
most of order (5^^ in probability by products of drift and Brownian increments in the overall correlation. □ 

In Figure[3]the observation times Tj, j = 0, . . . , 11 of F for our Example [T| from the last section are 
plotted against the observation times ti,i = 0, . . . , 11 of X. The dashed lines intersect for synchronous 
observation times to = ro,<3 = ts and tio = tiq on the diagonal of the square in Figure [3] A similar 
visualization of the realized covolatility estimator for synchronous and equidistant data would yield coex- 
tensive squares around the diagonal, over which multiplied increments are summed up. Refresh times are 
(in general) not equidistant but provide a synchronous realized covolatility estimator as an approximation. 
The Hayashi-Yoshida estimator (|9]l is the sum of products of increments with overlapping observation time 
instants. The relation to the synchronous approximation is that we have next-tick interpolations and 
previous-tick interpolations to the times Ti,i = 0, . . . , 8 and take increments from previous-tick interpo- 
lated values to next-tick interpolated values. The time instants of are visualized for our example in 
Figure [3] The previous- and next-tick interpolations are illustrated in Figure |4] The products of time in- 
stants leading to the error A^ are illustrated in the same picture by the grey rectangles. As can be seen 
for the example in Figure |4] A^ is the sum of the errors by the ith next-tick interpolation multiplied with 
the increments of the other process over [min {k, Xi), Tj] and the sum of the errors of the ith previous-tick 
interpolation multiplied with the increments of the other process over [Tj_ i , T^] . The sum of the increments 
over the squares in Figurejs] for our example, and the grey rectangles in Figure|4j A^ for our example, 
is the Hayashi-Yoshida estimator evaluated at the end of the last section. 
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Figure 4: Illustration of the next- and previous-tick inteipolated values and the error due to non-synchronicity for our example. 



Definition 2 (quadratic (co-)variations of time). For any G IN let Tj^^\ i — 0, . . . , be the times 
from the partition of [0, T] defined in ( fTT) above and gl^-* , 71^' , l^^^ , a|^' the corresponding observation 
times designated by Algorithm [T| from the estimator (|9]|. T/N is the mean of the time instants AT^^-* = 
j^(^) — T^^^ , i = I, . . . , N. Define the following sequences of functions 

G"W = f E {^Tr)\ (14a) 

T<"'<t 

+at;V"^ {rr - CO + Ai;^;^ {rr - , (i4b) 



T}"'<t 



^"(0 = ^ E (^^"^ - CO (^^^^ - ^^"0 + (^^"^ - (^f ^ - Tr) , (140 

for t E [0, T] that we call sequences of quadratic (co-)variations of times. 

A stable central limit theorem for the estimation error is deduced on the assumption that the sequences 
defined by ( |14a| i, ( |14b| l and ( |14c[ ) converge pointwise and the sequences of difference quotients uniformly: 

Assumption 3 (asymptotic quadratic (co-)variation of times). Assume that for the sequences of sam- 
pling schemes and the times gi^\ ^i^^^^^"^ '^^^ sequences of quadratic (co-) variations 
of times (t), (t), (t) defined in Definition^the following holds true: 

(i) G^(t) ^ G{t) , F^{t) F{t) , H^{t) H{t) as N ^ oo, where G{t), F{t), H{t) are 
continuously dijferentiable functions on [0, T]. 



13 



(ii) For any null sequence {hjy), hjy = O (^N ^) 



= 0{N 
















~F^(t) 




h-N 






h-N) 







G'{t) (15a) 
F'{t) (15b) 
H'{t) (15c) 



uniformly on [0,T] as N ^ oo. 



Assumption |3] is necessary to ensure that the sequence of variances of the estimator (|9]) converges as 
71, m — > oo. The derivative of the asymptotic quadratic variation of refresh times ( |15a| i will appear in the 
asymptotic variance of the discretization error , since refresh times are (in general) not equidistant. For 
^rj.{N) ^ rj^ij^ j^j. aiH e {0, . . . , N}, G'{t) = l[o,T] holds true. 

The uniform convergence of the difference quotients defined by ( |15b| i and ( |15c| i are necessary to ensure 
that the sequence of variances of converges as — oo. The assumptions imposed by ( |15a| i-( [T5cj i are 
weaker than assuming convergence of the joint sampling design of (7""^'" , T^'™) and are not very restric- 
tive. They hold true whenever the sequences of sampling schemes tend to a certain state of asynchronicity 
or have a uniform behaviour of non-synchronicity in the limit as rt, m — oo. For homogeneous sampling 
schemes these (co-)variations of time converge to linear limiting functions. 

The sequence of functions F^ describe an interaction of interpolation steps between the two processes. 
In contrast, Fl^ is defined to measure an impact of the in general non-zero correlations of next-tick and 
previous-tick interpolations to the same refresh time T,, for each process separately. 



Example: 

Consider the synchronous equidistant sampling schemes with N = n ^ m and = r|"^ = i/n,i 



0, . . . ,n. The left-hand side of Figure [s] shows the quadratic (co)variations of time G^ , and for 
N = 30000. F^ and i/^ are identically zero since there are no asynchronous observations and because 
T^^'' — i/N, = t/"' , < i < n, interpolation steps are redundant and equals zero. The function 
G^ is a step function that will tend to the identity on [0, T] as — > cx). 

Next, we consider a situation which originates from the complete synchronous equidistant one by shifting 
one time-scale half a time instant 1/2N. Then we have completely non-synchronous sampling schemes 
and we will call this situation intermeshed sampling. In this case the synchronous approximation is still 
equidistant with instants and, hence, G is the identity function. F and H are linear limiting functions 
with slope 1 and 1/4, respectively. Interpolations are carried out for all 1 < i < for the same process for 
which its first observation takes place after the first observation of the other process. All interpolation steps 
equal 1/2N and thus H' = 1/A follows. Since for H interpolated time instants 1/2N are multiplied with 
refresh time instants in both addends due to the specific structure, F equals the identity on [0,T]. The 
functions G^ , F^ , for intermeshed sampling are illustrated in Figurejsjon the right-hand side. 

In the next section we will show that for an important special case, independent homogeneous Poisson 
sampling, ( |15a| i-( [T5c] i are fulfilled when replacing deterministic convergence by convergence in probabil- 
ity. Furthermore, the stochastic limits G'{t), F'{t), H'{t) are calculated explicitly and are again constant 
on [0, T]. For data applications one can calculate easily empirical versions G'^ P'n 7n{t)jHn m(^) of 
G'{t),F'(t),H'{t) and use those as estimators for ( |15a| i-( [T5cl i. 

The key result of this section is the following Theorem[2] The detailed proof is postponed to the Appendix 
This result gives insight into the asymptotic distribution of the Hayashi-Yoshida estimator. It improves 



on the asymptotic normality result in ||T41 . since the weak convergence is stable in the setting where we 
allow for random correlation, drift and volatility processes. The representation of the asymptotic variance 
using ( |15a| i-( [T5cl ) differs from that in ifTSl . where a similar stable convergence result is established, by 
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Figure 5: Quadratic (co-)variations of time for synchronous equidistant (left) and intermeshed (right) sampling. 



the decomposition of the estimation error in ( [T2] i and ([T3| and the notion of (co-) variations of times. The 
latter provide helpful tools to describe the stylized facts and features of non-synchronous data and build 
the ground work for combined approaches for widespread generalizations and extensions of the underlying 
model. A very important one is the generalized multiscale estimator in fS] when market microstructure 
noise effects are taken into account. 

Theorem 2. The estimation error of (|9]) converges on the Assumptions^ ^and^stably in law to a centred, 
mixed Gaussian distribution: 

Vn (Xg, - XiJ (r^, - y^J - [X , F]^ j 4 N (0 , + VA^) , (16) 
with the asymptotic variance 

vu^+VA^ ^Tj^G'{t){a^arf{p1 + l) dt + T (f' {t){a^ ajf + 2H' {t){pta^ ajf) dt 
where the two addends come from the asymptotic variances of and A^, respectively. 



5. Independent Poisson sampling 



In this section, we consider the model in which the sequences of observation times are supposed to be 
realizations of two homogeneous Poisson processes that are mutually independent and independent of the 
processes X and Y. 

Thereto, let and m'-'^^t) be sequences of two independent homogeneous Poisson processes with 

parameters Tn/9i and TnjQi (n e M), such that the waiting times between jumps of n'^"^ and m*^"^ are 



exponentially distributed with expectations E 



At, 



(n) 



Qxln and E 



At, 



(n) 



Q2ln ,i e IN,j e IN. 



Thus, n^"^ (T) and to'^"' (T) correspond to the sequences giving the numbers of observation times of X and 



15 



Y in the time span [0, T]. The increments of the sampHng times of the closest synchronous approximation 
([TTJ are maxima of the exponentially distributed waiting times and we obtain: 



at; 



(n) 



F{t) = 1 - exp 



tn 



exp 



tn 

"e~2 



exp 



-t7l 



1 



, k e 



Denote N{T)^"^ = max^vgiN {X^a^o ^-^i"^ — -^}- focus on the characteristics of the sampling 
schemes affecting the asymptotics of the synchronized realized covolatility estimator (|9]l. In particular 
our interest is in the quadratic (co-)variations of times defined in Definition |2] 



Proposition 5.1. In the independent homogeneous Poisson model for sampling schemes, it holds true that 
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Proof. Poisson processes are Markovian and the exponential distribution of the increments between arrival 



times is memoryless. Wald's identity ensures that E 



l^k=0 ^-^t 



(n) 



E 



N{T) 



(n) 



E 



at: 



(") 



For 



the proofs of these attributes and further information on properties of mutually independent homogeneous 
Poisson processes we refer interested readers to Q- 

First of all we ascertain that ^ -rj"W(i,j) € {1, . . . , n(")(T)} x {1, . . . , m(") (T)} almost surely. 
For an arbitrarily fixed i, the expected values of next-tick, previous-tick and refresh time instants yield 
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The conditional expectations given that the ith refresh time T'^ 
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7j' ' is an arrival time of m^^"^ yield 
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r^(ri) _ ,(n) 
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since the latter previous-tick interpolation is zero with probability 1 if T;^"^ 7^ 7^"^. Only for (7;^"^ - A,^"^) 
the conditional expectation differs from the unconditional and can be calculated by further conditioning 



E 



rp(n) _ ^(n),rpin) _ (n) 
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(n) rp(n) _ , (n) 



+ E 



min) _ \ {n),rp(n) _ in) rp{n) _ ,(n) 



) 



= {91+92- 



0i92 



+29,-, 



1+^2/^1+^2 ^1 + t72 

where the factor 26i in the second addend is simply the expectation of the waiting time for two jumps of 
n. Here, we have used some simplifying symmetry aspects, a rigorous proof using the density functions is 
obtained by calculation of 

2^1 ^2 
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The conditional expectations on T> = gl are deduced analogously. Since E 



j,{n) _ j{n) 



E 



rp(n) _ rp 



(n) 



i-1 



-E 



rp{n) ,{n) 



and the (conditional) expectations of the products occurring in ,F^, 

H'^ equal the products of (conditional) expectations thanks to the memorylessness of exponential distri- 
butions, the latter results suffice to apply the law of large numbers to the empirical (co-)variations of 
times. For the asymptotics of G^{T),F^{T) and {T), we conclude for the number of addends 
7V(T)("), that EiV(r)(") = {T/9)n + o{n) with 9 = 9, + 62 - {9i92)/{9i + 6*2) what follows from 

E7V(T)(")E AT^^^ =T + Op{n-^) and Var (iV(r)(")) = 0{n-^) since 



Var I J2 AT^"M =Var(7V(r)("))E (at^^"))^ 



k=0 



+ E 



N{T) 



in) 



Var 



The exact probability mass functions of the counting processes N{t)^"^ associated with the maxima of the 



waiting times At^"^ , Arj"'' have a quite compUcated form, so that we only give the last two results on the 
expectation and the variance that are necessary for the proof of the proposition. 
From the preceding conclusions, it follows that 
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Figure 6: Quadratic (Co-)variations of times for liomogeneous Poisson sampling. 



Inserting 9 we obtain formulae ( |17a| i-( [T7c] i. In the evaluation of we have also used the second moment 
of AT^"^ which can be calculated using the above given distribution function. □ 

Figure |6] depitcs the quadratic (co-)variations of times for simulated mutually independent homoge- 
neous Poisson processes. On the left-hand side both parameters have been set 6* = 1 for T = 1 and 
n — 30000. The stochastic limits are linear increasing functions on [0, 1] with slope 14/9, 10/9, 2/9 and 
1/4, respectively. On the right-hand side we see the (co-)variations of times for T = l,n = 30000, 9i = 
1, 6*2 = 0.5. Those tend in probability to hnear hmiting functions with slope 82/49, 44/49, 8/49 and 2/9, 
respectively. 

In the model of non-synchronously observed Ito processes X and Y which fulfill Assumption [T| and ob- 
servation times following an independent Poisson sampling scheme of the above given form, we derive the 
following stable central limit theorem as special case of Theorem |2] 

Corollary 5.2. The estimation error of the synchronized realized covolatility estimator (|9| converges on 
the Assumption |7] conditionally on the independent Poisson sampling scheme with < (?i < oo and 
< ^2 < oo stably in law to a centred mixed Gaussian distribution: 




) (y („, - y,(„,) - [X , y]^ 4 N (0 , vt) , (i8) 



with the asymptotic variance 



where the two addends come from the asymptotic variances of the discretization error of the closest 
synchronous approximation ( |12| l and the additional error due to interpolations ( |13| l, respectively, and 
9 — 9i -\- 92 
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Proof. It is a basic result in the theory of extreme values that for the supremum of n i. i. d. exponentially 
distributed waiting times AT^ with EAT^ = it holds true that sup^ (AT,) = Op (log (n) /n). We refer 
to ifTTl for a proof. In the setting of mutually independent homogeneous Poisson processes with param- 
eters Tn/Oi and Tn/^s, we conclude that sup^^^^^ n{tY^)} ^ Op (^log 7V(T)(")/iV(r)(")y Hence, 
Assumption |2] holds for the sampling design where the orders of (5;f , hold in probability. Then all 



findings in the proofs of Propositions A. 2 and A. 5 stay valid when we insert the (co-)variations of time 



deduced above in the limits of the variances. □ 

The stable convergence holds conditionally given the observation times, what means that endogenous 
observation times are not covered but Poisson sampling independent of the processes X and Y . 
The asymptotic variance of the mixed Gaussian limit is in line with the results by 1 14] and |T5 1. We remark 
that one has to pay attention to the proportionality to 9 in the rate N[T)^^'^ when comparing the asymptotic 
variance to the one in ifTSll . 

From an applied point of view, the model considered in this section could be criticized for its flaw that 
sampling schemes of two correlated processes are modeled to follow two independent processes and for 
time homogeneity. Both seems to be rather unrealistic in financial time series. However, independent and 
homogeneous Poisson sampling times designs constitute the most commonly used model in this research 
area (cf. |23| , [ 13] among others) because they are handy and allow for explicit calculations while the 
model is not too far away from the real world. 



6. Asymptotic variance estimation 

Finally, we state a consistent estimator for the asymptotic variance of the Hayashi-Yoshida estimator 
Q from Theorem [2] Since in |14| a central limit theorem for the case of deterministic correlation and 
volatility functions has been proved, the asymptotic variance is non-random in that setting. In a recent 
publication |15|, in that the authors also generalize the asymptotic distribution result to a stable central 
limit theorem in the setting of random volatility and correlation functions, a consistent estimation method 
for the asymptotic variance is provided using kernel estimates. Our estimator differs from this method and 
we incorporate only one time transformed histogram-type estimator. 

Proposition 6.1. Define the estimator 

N-l 

AVAR^y E - Xi,){Y,^ - Fa,) [(^g, - Xi^W^^ - Y^^) 

+ 2(X,^.,, - X;,,J(K,,,, - - iTh 

with 




being a histogram-based estimator for {pta^ aYY'G' {t)dt. The estimators for the increase of the 
quadratic covariation on bins are Hayashi-Yoshida estimators of the type 



HY 



It holds true that 



AVARhy^T G 



'it) {a^aj)\pl + l)dt + Tf [F'{t) {afajfdt + 2H'{i) {p.a^ajfdt) 
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on the Assumptions^^and^ Thus, we have on hand a consistent estimator for the asymptotic variance 
of the Hayashi-Yoshida estimator and the feasible stable central limit theorem 

^^'^Jt A*n(0,1). (19) 



AVARffy 

For more motivation and details on the construction of histogram estimators for which bins are cho- 
sen equispaced according to a transformed timescale associated with a certain monotonic function, as the 



asymptotic quadratic variation of refresh times here, we refer to [5]- Proposition 6. 1 is proved in Appendix 

m 



A. Proof of TheoremU 



A.l. Discretization error of the synchronous approximation 

Proposition A.l. On the Assumptions^ ^and ( |15a| l the discretization error of the closest synchronous 
approximation converges stably in law to a centred mixed Gaussian distribution: 



0, 



G'{t) (afaf) + 



(A.l) 



Proof. In the proofs superscripts of the sampling times are frequently omitted to increase the readability. 
First note that on Assumption[T] by Girsanov's theorem we may without loss of generality further suppose 
that pf^ = p^Y = identically since we have learned in Section |2] that stable convergence is commu- 
tative with measure change. Let Mt and Lt be the continuous martingales Lt — /p dW^ , Alt — 
/p (tY dWj where W-^ , are two standard Brownian motions with quadratic covariation [W^ , W^^ - 

X 



pso-f erf ds and denote L, = erf dW^ , AI, = erf dWj 



Proposition A.2. On the same Assumptions as in Proposition A.l the process T)^ defined by 



N 



for < t < T converges as N ^ oo stably in law: 



V. 



N St 




(A.2) 



where W is a Brownian motion independent of T and 



(A.3) 



Proof. We will prove this stable convergence of the process associated with the transformed discretization 
error by application of Jacod's stable limit Theorem [T| It is also possible to use the discrete-time version 
of this Theorem from Corollary |2.3| which we apply in the next subsection. 

Using the definition of the quadratic covariation process of martingales or integration by parts formula, 
we find an illustration of the discretization error by a sum of stochastic integrals and an asymptotically 
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negligible term: 



T<"'<t 

= ^ {L,M,-L,M,_j 
r<"><t 



where we denote T(t) := maxi {T^^^^ < t). 
Thus, we obtain 



Consider the centred continuous martingale 



since [L, M]^^^^ ™ [L, M]^^ = [L, Af], + 

pT pT 



J%(t) Ji{t) ' 



We calculate the corresponding quadratic variation process at time t: 



_ N 
t ~ T 

N 



E 



T'"><t 



[L, - L,^i)dMs + / {Ms - M,^i)dLs 



%{t) 



T<"'<t • 



' [Ls-L,^ifd[M]^+ I ' [Ms- M,^ifd[L] 



Ti- 



%{t) 



N 



E 



+ 2 ^ " [i - [Af - M,_i], d [Af , L], ) + 

+ 2 / * [L - U_^l[M - M,_,ld[M,L]\ + Opil) 
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r 

n(N) 



T("'<t 

Jo 



In this calculation we have used integration by parts and the change of variables Theorem for the integrals 
with quadratic covariation integrators that are of finite variation. The second last equality is an applica- 
tion of the mean value theorem (the volatility and the correlation processes are continuous and thus also 
bounded on compact sets) where the constants {(J^)i, (o'^)i and {p<T^cr^)i come from. The Riemann 
sum converges and with Definition |2] and Assumption |3] this yields the convergence in probability of the 
quadratic variation to G'{s){p1 + l){a^ cr^ Y — /J w-p^. The third equality above is proved in: 

Lemma A.3. It holds true that the approximation error terms 

( [ ' (Ls- L,_ifd [M - M,_i]^ - / " [L- L,_i], d [M - Af,_i] J (A.4a) 

V ( / ' (M, ~ M,;-i)2d \L - - / " [M - d[L- J (A.4b) 



( f \m,^ M,_^){L,- L,_^)d[M,L\^- I ' [Af-M,_i,i-L,_i]^d[M,L]J (A.4c) 



T.<"'<t 



T 



converge to zero in probability. 



Proof. The proofs for ( |A.4a| i and ( |A.4b| i are completely analogous and we restrict ourselves to prove it for 
( |A.4a| i. By Ito's formula 

[L, - = 2 I {Lr- L,-i)dLr + [L- L,_i]^ 

holds. The left-hand side of ( |A.4a| i equals 

= J2 ( 2 / " {L, ~ U_i)i[M ~ M,_i]^J dL, - 2 I ' (Ls - U_i)i[M ~ M,_il) dLA 



T.<"'<t 
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by application of the integration by parts formula in the way 

Zt. [M - M,_i]j,. ^ r Ztd [M - M,_i], 



with := 2(Ls — ii_i)(iLs for Tj-i <<< T to the addends. Therefore, we can write the left-hand 

side of ( |A.4a| i in the way m''-^^ + -^2^' with two centred continuous martingales A^^^\ A^2^' defined 
in the fashion of 0'^-* above and calculate the quadratic covariation processes at time t: 



M 



[N] 



4 E 



T^"><t 



.fd[L]\+Op{l) 



T,-i 



< 



2 f'^' 

4max sup (L^ — ij_i)^max sup [M — Mi_i]^ / d[L]^+Op{l 



)■ 



The first addend is up to a logarithmic factor Op{Spf) and hence A^2^'' = Cp(l) on Assumptionji] That 

M.^^^ = Op{l) is proved analogously. This implies that ( |A.4a| i is C'p(l). 

The strategy of the proof for ( |A.4c| i follows the same approach, starting with the equation 



(L,-L,_i)(Af,-M,_i) 



Ti-i 



{Mr ~ Kh^i)d{L ^ L,^l)r 



Ti-i 



and applying integration by parts as above with Zt = (is — Li^i)d{M — Mi„i)s + {Ms — 
M,_i)d{L - L,_i)s for T,_i < i < T,. 

We complete the proof of the convergence of the quadratic variation with the proof for ( |A.4d| i. Denote 

{pa^a^)^ = (o'^)2j • (cr^)^j • to distinguish between the values from the application of the mean 
value theorems to the two different addends. An upper bound of the left-hand side of ( |A.4d| l can be found 
by elementary algebra and the triangle inequality for the absolute value: 

f Y: iATO'({p^a^r^+{-a^)Ua^r^-({pr^_,a^^_^^^^^^ 



T}"'<t 

N 



2 /,2 



T'"><t 



Op{l) ■ 



□ 



The martingales (p^^'> can be written for every N as time-changed Brownian motions B 



{DDS,N) 



by the Dambis-Dubins-Schwarz theorem. The sequence of martingales ^'^^ or associated time- 
changed Dambis-Dubins-Schwarz Brownian motions converges weakly to a limiting Brownian motion by 
the asymptotic Knight-theorem. The limiting Brownian motion will be defined on an orthogonal extension 
of the original probability space. To obtain the stable convergence result, we apply Jacod's Theorem[T]and 
thus, we are left to verify conditions ( paj i and ( pbji . 

Consider the quadratic covariation process of and the reference martingale L 



E 



Ti 



Ti- 



{Ls - L,_,)d [M, L]^ + / {Ms - M,^,)d [L]^ + Op{l) 



Ti- 
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The term of smaller order than 1 in probability comes from the increment of the covariation process on 
['X(t), <]. As before, this equality holds true for all t, since for t < Ti the covariation is Op{l). Integration 
by parts yields: 



E 



T'"><t 



([Af,L]^^ - [M,L]^^_J(L,-L,_i) 



Ti 



[M,L]^d{Ls-L,^i) 



Ti- 



Ti 



It remains to show that this term converges to zero in probability. The term is centred and using Ito isometry 
we find the following upper bound for the second moment: 



E 



<2^E 
- J, 



max sup ([M, L], - [M, L]^^ ^ 



+ max sup {[L]^-[L]^ d[M - M,,,]^ 

je{l,...,JV} sg(T._i,Ti] ^ JTi.i 



Ti 



O {NS 



N) ■ 

The term is bounded by a constant times NSj^ since squared increments, cross products of increments 
and increments of the quadratic (co-)variations of L and M over time instants AT^^'' are bounded by 
At/^'' times a constant. To sums with products of time instants we can apply Holder's inequality with 
the supremum norm to obtain upper bounds. There are at most order S]^^ time instants AT^'^^ of order 
sup, AT;^^^ = Sn since J^t ^Tj;^^ < T and the time span T is fixed. 

Hence, [i,^^^^]^ = Vi G [0,T]. With the same strategy [M,0(^)]^ = Op(l) Vi G [0,r] can 

be shown. 

For every bounded J^t -martingale satisfying [i, L-'-] = the covariation 



N ^ 



Ti-i 



{Ljr - Lt,)d [M, L-^]\+ = 0,(1) 



converges to zero. The same holds true for every bounded J^t -martingale orthogonal to M. Applying 
Theorem[T] we deduce that Proposition A. 2 holds true. □ 



Proposition A. 1 is a direct consequence of the stronger result in Proposition A. 2 since for t = T the 



marginal distribution is a mixed normal distribution which is independent of T. The stable convergence 
assures that the convergence also holds under the original probability measure and non-zero drift terms 
with the same asymptotic law. □ 



A.2. Error due to non-synchronicity 

Proposition A.4. Let Assumptions [7] |2] and ( |15b[ )-( [T5c] l from Assumption |i] be satisfied. The error A, 
due to the lack of synchronic ity converges stably in law to a centred mixed Gaussian distribution: 



(A.5) 
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with asymptotic variance 

VA^ = r F'{t) {afajfdt + r 2H'{t) {ptafajfdt . (A.6) 

Jo JQ 

Proof. First, we write the ith increments occurring as factors in the addends of the estimator (|9]) as the 
sum of the next-tick interpolation at Ti, the increments AXxi — Xt^ ~ -^Ti_i and AF^; — Yrt — Yxi-i, 
respectively, and the previous-tick interpolation at Ti^i and multiply out the addends. 

N 

i=i 



N 



((Xg^ -XtJAFt. + (y^, -FtJAXt. + {Xt,_, -Xu)/XYt, + {Yt^_, - IaJAXt, 

+ (Xg, ~ XTMYT^-^ - + {Y^, - YTMXT^-^ - X^)) + 

The indicator functions in ( [T3| ) have been dropped since the corresponding addends are zero if the indicator 
functions were zero. Since at least one of the next-tick interpolation errors is zero and as well one of the 
previous-tick interpolation errors, too, two addends, namely the products of next-tick interpolation errors 
and the product of previous-tick interpolation errors, equal zero. Thus, the error due to asynchronicity can 
be written as the sum of the remaining six terms (where at least another three equal zero in each addend). 
We conclude, that the error can be expressed in the following way: 

N-l 
i=l 

{Xt,^, - XtMYt^ - Fa.+J + {Xt^ - Xi,^,){Yt^^, - Ft.)) + . 

In this equality an index shift has been applied to the partial sum of previous-tick interpolated errors multi- 
plied with AXxi and AFy. , respectively, leading to the structure that in the ith addend the factors contain 
next- and previous-tick interpolated errors to the same T^. The Op(l)-term emerges from end-effects when 
shifting the original sum. 

In the last illustration of A^ consecutive addends of the sum are uncorrelated in contrast to the non-shifted 
illustration. The reason is that, if without loss of generality ji — Ti holds, {Xg^ — XtJAFt-^ and {Xxi — 
X; ) Alj-.^j have in general a non-zero correlation whereas {Xg^ — Xy. ) Al^. and ~ Xi.)AYt- 

are uncorrelated. Furthermore, the fact that = Ti ^ A^+i = Ti assures that the addends in the last 
illustration of A^ are uncorrelated. Roughly speaking we capture correlation between subsequent addends 
of the outer sum and transfer it into additional correlation in the inner sum. 

As in the foregoing proof of Proposition | A.l[ it is sufficient to prove the stable convergence result for 
the zero-drift case. We denote, as before, the corresponding transformed processes Lt = /g dW^ and 

M, = J^aYdWr. 
Consider the sum 

= ^ AAf ^ ((Lg.-LTj(MT,-M,J + (M^.-MT.)(LT.-i/.) 

-"i + l — i + 1 — 

+(Lt. - Li,+J(MT.+, - MtJ + (Mt. - Ma,+J(Lt,+i - ^tJ) (A.7) 

for fixed < i < T. 



Proposition A.5. Assume the same conditions as in Proposition A.4 For fixed < t < T the transformed 
error due to non-synchronicity is the endpoint of a discrete, centred, square-integrable martingale with 
respect to the filtration J-i n :— J- (n). The process Af converges as N oo stably in law: 

Af 4 A = /* V^dW^ (A.8) 

-'0 
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where W-^ is a Bwwnian motion independent of T and 



Proof. The expectation of the absolute value of the sum is bounded for all t E [0,T] and A A. 
0, . . . , iV are J^i^N 



(A.9) 



N 



J"„(]v) -measurable. Since 



= E [Lg^ - Lt,] {Mt, - Mx,) + E [M^^ - Mt^ {Lt, - 

+ {Lt, - i/.^JE [AA/t.+J + (AfT. - AfA.+JE [ALt.+J = 

for the conditional expectation of the increments holds, Af' is the endpoint of a J^i ^v -martingale. 



The stable weak convergence to a limiting Brownian motion is proven with Corollary 2.3 to Jacod's Theo- 
rem [T] 

First, we verify the conditional Lindeberg condition that is implied by the stronger conditional Lyapunov 
condition. It is sufficient to proof the following: 

Lemma A.6. The sum of the conditional fourth moments of the martingale increments Af^ converges to 
zero in probability: 



E 



N 



T"2'<t 



= Op(l). 



Proof. Throughout the proof C denotes a generic constant that does not depend on N . We consider dif- 
ferent addends of the fourth conditional moments consecutively. The sum of conditional fourth moments 
incorporates addends of the following types: 



• fourth-order moments: 



2^2 



J2 E [{Lg^ - LtJ^] (Mt, - A/a J' 



T.^:v<t 



• second-order moments 

iV2 



y2 



T "'<t 



• third- and first-order moments 

7V2 



rp2 



For the partial sum with addends of the first type an application of the Burkholder-Davis-Gundy (BDG) 
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inequalities yields 

2^2 



J2 E[{Lg,-LTj']{MT,-M^J' 



T}"'<t 



<c 



J^2 



T"'<t 



Ti 



[Mt^-MxJ' 



N' 

<C-^ sup (af )2 J2 {MT^-M^y{g.,-T,f <Op{NS%)^Opil). 



rp2 



T"><t 



The last inequality can be deduced by the result that the convergence (A^/(3T)) ^^(AMtJ^ — ?> J^icrY^ds 
holds almost surely as iV — > cx) for the so-called realized quarticity (|3|) and that (g^ — T^) < Sn- Without 
the result about the convergence of the realized quarticity, the asymptotic order in probability can be de- 
rived by the convergence to zero of the expectation of the above sum and calculating the second moment 
that is bounded from above by a constant times N'^6lf. 

For the partial sum including addends that incorporate second-order moments we obtain an upper bound 
by application of the Cauchy-Schwarz inequality and the BDG inequalities: 

6E [(ig, - LtJ'(AA/%^J2] (^^^ _ L,,^J'(A% - A/a J' 



y2 



T"'<t 



2^2 



T)' <t 



T"'<t 



{affds 



Ti 



E 



{a^fds 



Ti 



Op(l) 



The stochastic order follows, since the term has the expectation 



C 



iV2 
2^2 



E 



{affds 



E 



{alYds 



<C^^6lE(/"(affd.)E 



2^2 



{(T^fds] E 




(af)^d.) El / {c^lfds 



Ti 



2\ 2 



where again the Cauchy-Schwarz and BDG inequalities have been applied. The variance is bounded from 
above by a constant times N*Sj^, what can be shown by a similar calculation where thanks to the fact 
that Ti — -fi Ai+i = Ti the addends are uncorrected and the variance of the sum equals the sum of 
variances. 

We treat the third type of addends occurring in the sum of conditional fourth moments in the same way. Ito 
isometry yields 



2^2 



T"'<t 



7V2 

2^2 



T"'<t 



Psufa^ds 
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This term has expectation 

Ar2 



^■2 



r<"?<t 



Ti 



< 



J'2 













E 









< CN^5% = o{l) , 

and an analogous calculation as before yields that the variance is of order N'^Sj^. 
Thereby, the sum converges to zero in probability. 

Next, we consider the sum of conditional variances of the increments of the discrete martingale. 
Lemma A.7. 



□ 



E 



4 /* F'{s) {a^aj)'ds + f 2H'{s) {psafa^fds . (A.IO) 
Jo Jo 



It holds true that 
Proof. 



E 



E i^^ry 



= ^ E ^[{Lg,-LT,f{MT,-Mx,f + {M^,-MT,f{LT,-Li,f + {LT,-Li,^,f{AMT,^,)' 



r(")< 



N 

(Misometry) 'J^ 



E E 



{a^fds 



(Lt. -ii.+J'E 



[a^yds 



Ti 



{LTi-Lif 



{a^rds 



+ 2(M-r, - MaJ(Lt, - ii,+i)E 
+ 2{Lt,-Li^){Mt,-Mi^^JE 



I psufa^ds 

jTi 
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N 

(LemmifXsl T 



(N)^ 
i + 1 i 



{affds 



+ / {affdsV. 

PsCyf ds E 











+ / ' [a^fdsV. 


r"\affds 









Psa^ajds 



Psa^a^ds 



-Op{l) 



N 



{affds I {alfds+ / {alYds I {at Yds 



Ti 



Ti 



E 

+ l^\atrds I {crlfds+l {alfds 



Ti 



+1 

Ti 



Ti 



Ti 



N 



Pscrfa^sdsl pscrfa^ds + i I psafaldsl psCrf a; ds \ ^ Op{l) 



T.("'<t 



+ (T, - A,+i)AT,+i) + (pT.CT^fT|:j2(2(r, - k+i){g, - T,) + 2{T, - A,+i)(7, - r,)))+Op(l) 



E 

T.<"'<t 



F{Ti+i)-F{Ti) xY'.2;.rp H{Ti+i) - H{Ti) „x y -.2 A,rp 



i+1 



i+1 



F'(.)(afar)^d.s+ / 2i7'(,s)(p,afar)^ds 



In the last step we have involved Definition |2] The Riemann sum converges on the Assumption |3] (in par- 
ticular \\5b\ and ( |15c| l) in probability as — > cx) to the expression /* vj( ds with given in Proposition 
[A3] 

The detailed proofs of the approximations are postponed in the following two lemmas. 
Lemma A.8. On the assumptions as before, the following equations hold true: 









l^aJfdsjE 


r\afrds 


= Op(l) , 




JTi 





iV 



E 

T^S<t' 









l^aJYds (e 


' Tia^fds 






-JTi 





N 



E [i^T.-Li 



T'^'^Ht 









Uif- 


^\cjffds\E 


' Tiajrds 






Jt^ 



N 
f 









j\<jffds {e 




- f\alfds^ 









N 



T^"^<t 









/J' i^IYds^ E 


£'\<^ffds 
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N 



{erf Yds 



{afYds\=0,(l) 



T<"'<t 



, Jli+1 J JTi 



T 



JTi JTi J 



N ^ 



ds 





■ /-Si 











TV 



(N)^^ V + 1 / \-JTi 



T^S<t 



ds 



N 
f 



T}7'<t 



Proof. We restrict ourselves to the proof of the first two equalities, since all other terms can shown to con- 
verge to zero in probability in an analogous way. The left-hand side of the first equahty has an expectation 
equal to zero which can be concluded directly by Ito isometry: 



E 



N 







j\alYds\E 


f\affds 







. 



In order to derive the stochastic order of the term, consider the second moment: 



E 



LV 

p2 



N 
T 



E ({MT^-M,^r~£\arrds]E\^£ 



E E 



{a^fds 



T"'<t 



X E 



{affds 



o{l) 



where the asymptotic order is deduced by Ito isometry and the BDG inequalities. Since the error induced by 
this term in the approximation of the conditional variance before is centred and has a variance converging 
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to zero as — > oo, the error is asymptotically negligible. 

In the second equality we consider the error when the expected increment of the quadratic variation of X 
over the next-tick interpolated time interval is substituted by the integral itself. We proceed as before for 
the first approximation. Since 



E 



N 
T 



Ti 



T'"'<t 



{affds 



Ti 



Ti 



= 



and 




E / '{-"^fdsU 



T "'<t 



{afYds 



{affds 



-oil), 



the approximation error is asymptotically negligible. The fact that ji ^ Ti ^ A^+i = Ti has been used 
that guarantees that the addends of the sum are uncorrected. □ 



Lemma A. 8 has been applied in the second and third equality in the evaluation of the sum of conditional 
variances and the proof of Lemma A. 7 is completed by the following 

Lemma A.9. On the same assumptions as before, the following equation holds true 



N 
T 



E 



T'2-,'<t 



Ti 



Ti 



and analogously the errors in the five other addends converge to zero in probability when replacing the 
product of increments of quadratic (co-)variations by the values o/py. , cr^ , a^, multiplied with the corre- 
sponding times increments. 

Proof. We prove the equality explicitly given in the lemma. The five remaining terms can be handled by the 
same strategy. By an application of the mean value theorem, elementary algebra and the triangle inequaUty 
for the absolute values, we deduce that 



N 
T 



N 



T 



J2 ( f\^^fds l ^\aj fds ~ {.alal:)\T, - \){g, ~ T,) 



T"'<t 



N 



T'"'<t 



< 



E ff -p 



''f-i'^Tf\+ sup \{aJf-ialf\]m-X,)ig,-T,) 



holds on Assumption[T] 



□ 
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□ 



To prove the stability of the convergence in Proposition ] A.5[ we show in the following that the discrete 
covariations of Af^ with the J^-generating underlying martingales L and M converge to zero in probability. 



Lemma A.IO. 



n"'<t 



Proof. Both relations are proven similarly and we restrict ourselves to the proof of the first one. The 
left-hand side equals 



y (N) 



(A/t. - ) + E 



{LTi — Li, 



Ps<Js 



Ti 



=: r 



r is centred and using Ito isometry the variance is shown to converge to zero: 



TV 



E 



{alYds 



+ E 



Psafa^ds 



E 









+E 


r 











+2E 



-2E 



Ti 



Y Y J 



Ti 



E 



E 



{at yds 











j +E 


r i'^Yfds 


^E 


r"\atfds 








JTi 



H \ 2 



E 



E 



Ps(^s ds 



PsCTs <^s ds 



UT, 



< CNS% ^ o(l) 



Once more we can conclude that the addends are uncorrected since =Ti^ = Ti and gi = Ti ^ 
Zi+i — Ti, respectively. □ 

Finally, we prove that the discrete covariation of our considered martingale with every bounded Tt- 
martingale that is orthogonal to Lt or Alt, converges to zero in probability. Hence, this lemma wiU complete 
the proof of Proposition] A.5| 
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Lemma A.ll. Assume that and are bounded J^t-martingales, with [L, L^] — and [M, Af ^] — 
0, respectively. It holds true that 



i+l i:^ 



T":'<t 



J-'rp{N} 



0, 



Proof. As in the preceding lemma, we only prove the first part of the result. The left-hand side of the first 
equation equals 



J2 E [aL^^^^ ((Lg, - LtJ(Mt, - Ma J + (M^, - Mt.)(Lt. - LiJ 



<t 



a 



T<"'<t 



d[M,L^]^ 



— V E 

T<"'<t 



d [M, L-L] 



+ (A% - Ma,+ JE 



Ti 



d[L,L^]^ 



This term is centred and the has the variance 

-J-l 1 



T E E 



T.<")<t 



d [M, 



E 



-E 



7. X 2 



X E 



TV 

Ito isometxy T 



7,4 



d[Af,L^]^ 



Ti 



E 



d[Af,i^]^ 



d[Af,L^]^ 



E IE 



+ E 



■E 



d [M, 
d [M, L-L] 







E 






Jii+i 



{a^rds 



E 



E 



(a^fds 
d [M, 



E 



Ti 



d [A/, L-L] 



Thus, the covariations converge to zero in probabiUty. 



The Lemma completes the proof of Proposition A.5 



□ 
□ 



The mixed normal limit in Proposition A. 4 can be obtained as the marginal distribution of in 



t = T. 



□ 



Proposition A. 4 for the error of the approximation by the discretization error of the closest synchronous 
approximation ^3 1 and the stable limit theorem for this synchronous discretization error ([12]) given in 
Proposition ] A.l suffice to imply Theorem|2] That is because the multivariate stable convergence Theorem 
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|2.4| applies to the vector of the two uncorrelated terms and since the covariations converge to zero , the sta- 
ble convergence to the mixed Gaussian limit with the sum of the two asymptotic variances is concluded. □ 



B. Proof of Proposition |6ll] 



The proof will be divided into three parts in that the sum of squared products, products of consecutive 
increments and the histogram estimator are considered, respectively. Denote ~ Xg. — Xtj , Xj' — 
XTj_i — Xi- and Xj = Xtj — ^Tj^i , j — 1, ■ ■ ■ , N. In the first step it is proved that 



N-1 



X, 



7j 



,2 p 



/ G'(i) (af af)' (2p? + l)dt + T fp^t) {af ajf dt 
Jo Jo 



N-1 



N-1 



ny: {x+ +x^' +xry {y; +Yf +Yry^Nj2mriYf +Y-r+{Y;nx^ +x-r 

+{xvf{Yfr+{Y-rixfr+{x^Yf)')+o,{i) 

All centred addends have a variance tending to zero as oo and converge to zero in probability. 

The sum of the first four adde nds ti mes the factor N/T has been proved to converge in probability to 
F' {t){a^ Y dt in Lemma A. 7 where this term has appeared in the sequence of conditional variances 
of the error due to non-synchronicity. 

Hence, it remains to prove that N Y!^^^ {XfYf f 4 T {2p1 + l){a^ aJYC {t)dt. For this purpose 

recall the notation from the proof of Proposition A.l With it = jl^fdW^ , Mt = jlcf'^dWj, 
Li = , Mi — , we can write the term 

N-1 N-1 I T, 

^ V ((L-L,_i)t,(A/-M,_i)tJ' = A^ V 2 / (L-L,_i)t(M-M,_i)2d(L-L,_i)t 

j = l ^=l V 



where we have applied Ito's formula. The sum of the first two addends converges to zero in probability 
since it is centred and the variance converges to zero. Since 

'(L-L,_i)t(M-M,_i)td[M-Af,_i,i-L,_i], = / ' {L - U_i)t{M - M,_i)td[M,L\ , 



the sum of the third addends has been considered in the proof of Proposition A. 2 as part of the quadratic 
variation of the discretization error of the closest synchronous approximation and converges in probability 
to 2 T G' {t){pta^ Y dt. The remaining sum of the fourth addends is also similar to the other part of 
the quadratic variation in the proof of Proposition] A.2| An analogous approximation and integration by 
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parts yields 



[M-M,_,]^d[M]^ + [ ' [L-L,_i],d[i],+Op(l) 



d{[L - L,^,]^[M - M,^i],) + Op{l) 



and the convergence of the above given term to T {a:^ <tY G' {t)dt. 
In the second part of the proof we are concerned with the term 



Af-l 



2N {X+ + Xf + X-){X+^, + XjVi + XJ^,){Y+ + Y^^+ Y-)iY+, + Y^^, + Yf^,) 



N-l 



= 2^ E {xfYfX^^^fi, + X+YfX-^,Yf^, + Y+X^Y^^Xf^^ + o,(l) . 

3 = 1 

The sum incorporating all centred addends converges to zero in probability. The last two addends capture 
the only dependence between consecutive addends in the error due to non-synchronicity ( [T3| ), namely when 
next-tick interpolations and previous -t ick in terpolations at the same Ti,i = 1, . . . ,N are included. Those 
have appeared in the proof of Lemma A.7 and have been proved to converge to T J^^ 2H' {t){a^ crY)^dt 

in probability. That 2NJ2{XfYfX^^^^Yf_^^) 4 2 /J G'{t){ptaf a^^dt follows with the methodology 
from |5 | and Lemma 1 from [23J using the concept of a time-change in the asymptotic quadratic variation 
of refresh times such that J^ii^'^i ~ T/N)"^ — o{N^^) holds true. Using the mean value theorem and 
(AT-jf - ATjATj+i = ATj{ATj - T/N) + ATj{T/N - AT,+i) together with the Cauchy-Schwarz 
inequality 



N 



N-l 

E 



I (at, (at, 



T 

N 





JV-1 




E(AT.)\ 







N-l , 

^ ( AT, 



T 
N 



yields the result. 

The Hayashi-Yoshida estimators on the bins in the histogram-based estimator (fT9]l fulfill 



A\Xyt"r ^ [ ' ptafaJdt + oJKfN-'l' 



so that the estimation error of the sum is of order K^J^'^N in probability and for K 



OO, N ^ 



LATjv — ^ consistency holds and we conclude consistency of the estimator of the asymptotic 
variance. □ 
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